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Standard Big Bang nucleosynthesis at baryonic density as inferred by WMAP implies a primordial 
^Li abundance factor of two to three larger than that inferred by observations of low-metallicity 
halo stars. Recent observations of ®Li in halo stars suggest a pre-galactic origin of this isotope, 
which is exceedingly difficult to explain by putative high redshift cosmic ray populations. We 
investigate if one or both of these lithium problems may be solved by late decaying relic particles in 
the Constrained Minimal Supersymmetric Standard Model (CMSSM) coupled to gravity. Assuming 
that the gravitino is the lightest supersymmetric particle (LSP) we find that in large parts of the 
CMSSM parameter space both of these problems can be potentially solved. In particular, we find 
solutions where both lithium problems may be solved simultaneously. These entail the hadronic 
decay of relic C>(lTeV) staus into C9(100GeV) gravitinos at ^ 1000 sec after the Big Bang, as 
proposed by one of us before [Ij]. Moreover, the produced gravitinos naturally contribute a large 
fraction, or all, to the by WMAP required dark matter density. A prediction of this model is the 
dark matter to be lukewarm with free-streaming lengths of the order of a 4 keV early freezing-out 
relic particle. Such lukewarm gravitino dark matter may have beneficial effects to the formation of 
galactic structures and may be constrained by the reionisation history of the Universe. The current 
version of the paper presents results for both cases (a) when catalytic nuclear reations are included 
and (b) when they are neglected. 


I. INTRODUCTION 

Determining the nature of the cosmological dark matter is one of the main outstanding challenges in modern 
cosmology. Compelling candidates for the cold dark matter may be in the form of new neutral stable particles 
arising in phenomenolo^ally successful supersymmetric extensions of the Standard Model of particle physics. These 
include the neutralino 0, a mixture of the supersymmetric partners to the neutral Hig^, SU(2) W, and U(l) 
B bosons, the axion and its supersymmetric partner the axino Q and the gravitino [JQ (the supersymmetric 
partner of the graviton). Whereas the bulk of studies has been devoted to the neutralino, the case of gravitino dark 
matter has recently attracted much attention. Studies of gravitino dark matter have been performed either in the 
context of smersymmetry breaking in a hidden sector being communicated to the visible sector by gravitational 
interactions [1,0) iinaini or by gauge interactions [H| . Gravitinos may be produced in the early Universe by (at 
least) two mechanism: (i) scattering of thermal radiation at the highest temperatures of the early Universe, hereafter 
referred to as thermal production (TP) and (ii) freezeout and decay of the next-to-lightest supersymmetric particle 
(NLSP) to gravitinos, hereafter referred to as non-thermal production (NTP). Whereas the gravitino yield in TP 
depends on the reheating temperature Tr, the gravitino abundance due to NTP may in principle come by itself, and 
independently of Tr [T^I , very close to the range inferred for dark matter density from cosmological observations. In 
the case of the gravitino LSP, NLSPs typically decay during or after the epoch of Big Ban g nu cleosynthesis (BBN), 
unless the gravitinos are rather light, thereby potentially disrupting light-element yields [a, [ij. This has been often 
taken to disfavor for heavy gravitino dark matter. However, though stringently constrained by BBN, recent careful 
study shows that much viable gravitino dark matter parameter space remains i 0,113, Eli- 

The epoch of BBN has long been known to synthesize the bulk of the '^He and D, as well as good fractions of the 
^Li and ^He in the presently observed Universe. Paramount to the realization of this fact was also the discovery of 
the ^Li-” Spite” plateau in 1982 [13], in particular, the observation of constant ^Li/H abundances in low-metallicity 
Pop II stars over a wide range in metallicity. This indicated a pre-galactic origin of ^Li, as other sources (i. e. g alactic 
cosmic ray nucleosynthesis) predicted a rise of ^Li with metallicity. Current observational estimates [13, [iJ, [l3| of 
^Li/H range between 1.10 x 10“^° and 2.34 x 10“^°, with differences mostly depending on which effective stellar 
temperature calibration for the Pop II stars is used. With the accurate estimate of the baryonic density by WMAP, 
i.e. ~ 0.022 — 0.023, it was possible to predict the primordial ^Li/H abundance of ^Li/H ss 3.82 — 4.9 x 10“^° [l3| 
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within the framework of a standard BBN (SBBN) scenario. It is apparent that this predicted abundance is a factor 
of 2 — 3 larger than that observed. 

It is conceivable that ^Li in the atmospheres of Pop II stars has been transported down beyond the base of the 
convective zone of the stars, and thereby depleted by nuclear burning (i.e. '^Li(p, a)"‘He). Though standard stellar 
models may not account for this ^Li depletion within the near-turnoff, low-metallicity stars of the Spite plateau, 
effects not included in standard models, such as rotation, atomic diffusion, or gravity waves could potentially change 
this conclusion. Nevertheless, rotationally induced depletion of ^Li by a factor of 2 — 3 predicts a spread in stellar ^Li 
abundances [2^, not observed by any group (see, however. Ref. [11). This conclusion could be potentially changed 
when either internal stellar gravity waves [22l| or magnetic fields are considered in conjunction with stellar rotation. 
Atomic diffusion, a process required to understand the structure of the Sun by helioseismology, predicts a slo pe in 
the plateau (as a function of stellar temperature) [2^ which is not observed in the data (see, however, also Ref. [2^). 
Only when atomic diffusion is coupled with an ad hoc and fine-tuned weak turbulence in the radiative zone of the 
star, may a depletion of the required factor of 2 — 3 result [2^. However, even in this case a dispersion of the ^Li data 
is likely to emerge. It is thus not impossible that depletion of ^Li on the Spite plateau has occurred. Nevertheless, 
^Li depletion in low-metallicity stars is a well- and long- studied possibility, and even after about 20 years of efforts 
no consistent and well-motivated scenario including a factor of (^2 — 3 depletion has yet emerged. 

Standard BBN leads to the synthesis of a ®Li abundance of ®Li/H at the level of —10“^^, orders of magnitudes 
below what is observable by current technology. Observed ®Li abundances in the Sun, galactic disc-, and halo- stars 
are thus traditionally not believed to be of primordial origin, but rather due to galactic cosmic ray nucleosynthesis via 
supernovae produced energetic p, a oi CNO inducing spallation p + CNO LiBeB or fusion a -I- a ^ Li reactions 
in the gas of the interstellar medium. Here ®Li is produced along with ^Li,®Be, ^°B, and ^^B. Typical production 
ratios of ®Li /®Be ~ 5 — 10 in galactic cosmic ray nucleosynthesis are consistent with those observed in the Sun, but 
not with those ®Li/®Be 40 — 80 ratios observed in low-metallicity \Z] ~ —2 Pop H halo stars. It is thus clear that, 
at the very least, the composition of galactic cosmic rays at higher redshift has to be strongly modified in order to 
account for the observations. Nevertheless, a ®Li/H abundance ratio of ~ 10“^^ at [Z] ~ —2 was concluded to be only 
with difficulty, or not at all, synthesized by galactic cosmic rays, an argument based on the energetics of supernovae 
generated cosmic rays [^ . 

It had therefore been speculated that the ®Li abundance in Pop H stars may originate from the early Universe, via 
the electromagnetic decay of a relic particle, such as the gravitino, inducing the non-thermal nuclear reaction sequence 
of ^He( 7 ,p)^H photodisintegration with the resulting energetic ^H (and ^He) further fusing on "^He to form ®Li 
Here it was found that the synthesis of ®Li was efficient enough to produce the Pop H abundance without disturbing 
the other light isotopes or the Planck spectrum of the cosmic microwave background (CMBR) at the observable level. 
With the advent of the first fully coupled calculations of thermal nuclear reactions and cascade nucleosynthesis during 
BBN [I|, further solutions for the synthesis of the ®Li abundance in the early Universe were found. They were 
based on either the residual hadronic annihilation of a population of dark matter (e.g. neutralinos) during and 
towards the end of BBN, or the hadronic decay of a relic, long-lived particle population around 10^ sec after the 
Big Bang [I| (e.g. gravitinos). In hadronic decays or annihilations the initial photodisintegration reaction in the 
non-thermal nuclear reaction sequence given above is replaced by ‘^He(iV,p)^H spallation by energetic nucleons N 
produced during the decay. 

An intriguing further consequence of the hadronic decay of a relic particle at 10^ sec is the prediction of a significant 
^Li abundance reduction concomitant with the ®Li production. This occurs due to the thermal nuclear reaction 
sequence ^Be(n,p)'^Li and ^Li(p, a)'*He induced by the excess neutrons due to the decay. At the same time the D is 
increased due to p(n, yjD , but not as much as to exceed a conservative observational upper limit of D/H < 5.3 x 10“® 
on this isotope. It was found [l| that a relic particle abundance of 0,xh^Bh Ri I — 5 x 10“'^ (depending on the 
mass of the relic), where Bh is the hadronic branching ratio of the relic A, was sufficient to explain qualitatively 
and quantitatively both, the low observed ^Li abundance and the high observed ®Li abundance. In [I| it was further 
speculated that among other possibilities, the relic could be a supersymmetric stau of ~ I TeV mass decaying into a 
LSP gravitino of mass 50 GeV. A clear prediction of the proposed scenario is the existence of a ®Li plateau, analogous 
to the ^Li-Spite plateau, i.e. constant ®Li /H in low-metallicity stars. 

Earlier suggestions of a primordial solution of the ^Li discrepancy Q invoked an electromagnetic decay of a relic 
particle around 2 x 10®sec, thereby photodisintegrating ^Be. Synthesis of ®Li during the same process was not 
considered. It was subsequently shown that such a scenario may not work, as either the observational upper limit 
on the primordial ^He/D ratio is surpassed due to concomitant ^He photodisintegration, or a reasonable lower 
limit on the primordial D /H > 2.2 x 10“® due to D photodisintegration is violated. 

Over the last months the number of p reliminary detections of ®Li/H in Pop H stars has multiplied by a large factor. 
There are now around ten claimed detections of the ®Li/^Li isotope ratio, with all ratios falling in the range 

between ®Li/^Li « 0.03 — 0.07 (and with average ®Li/’^Li« 0.042), independent of stellar metallicity falling in the range 
\Z] ~ —2.75 and \Z] ~ —1.2. There exist also around ten upper limits, with all of the stars, nevertheless, consistent 
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with ®Li/^Li in their atmospheres on the level of ^0.01 [s^. It is thus intriguing that the ®Li data indeed shows 
a plateau within a large metallicity range. From a galactic cosmic ray point of view the high ®Li/^Li abundances 
as reported in the lowest [Z] ~ —2.75 metallicity star LP815-43 is particularly difficult to explain. Following the 
announcement of results of these difficult observations with the VLT telescope, preliminary detections of more ®Li/^Li 
ratios by the Subaru telescope and Ref. [s^ have been claimed, with, in particular, one star of metallicity [Z] ^ —3.25 
seemingly showing again a similar ®Li/^Li ratio [1^. All this data points to the existence of a ®Li plateau (see 
below, however). Nevertheless, it needs to be stressed that each individual detection of ®Li/^Li due to these difficult 
observations is only at the about 2-4 sigma level. In this sense the observations have to be taken as preliminary. 

Though the observations seem to indicate a ®Li plateau over a wide range in metallicity (in particular also, when 
disk stars at metallicity [Z] ^ —0.6 are included into the sample), a plateau may not necessarily exist if one is to 
believe the claimed [4^ metallicity dependence of ®Li destruction on the pre-main-sequence (PMS) of the observed 
stars. Predicting ®Li destruction during this phase is less certain [4^ than du ring the stellar main sequence. In fact, 
the prediction are not comparing well to observations and are therefore suspect ^|. The predicted ®Li destruction 
would imply a rise of ®Li with metallicity for stars with [Z]^ —2. This indeed is theoretically favorable as most groups 
also observe a rise of ^Li with metallicity on the Spite-(quasi)-plateau, usually attributed to cosmic ray production 
of ^Li. As ®Li is produced by the same process a consistent picture would emerge. In the absence of PMS ®Li 
destruction either a fine-tuned stellar main-sequence ®Li destruction would have to be invoked [4^, or the conclusion 
of ^Li rising with metallicity would be erroneous. Here the latter possibility seems likely as the existence of a slope 
in the Spite-plateau is not confirmed by all groups. In contrast, though further observations are required, stars with 
metallicities [Z] ^ — 2 seem currently, in any case, to be consistent with a low-metallicity ®Li plateau, even when 
stellar PMS effects are included. 

The new ®Li data has already prompted the first attempts to be explained in terms of cosmic ray nucleosynthesis. 
In Ref. cosmic ray nucleosynthesis with energetic a’s generated at the shocks resulting in merger events during 
the formation of the Milky Way were proposed to possibly generate the ®Li abundance of ~ 10“^^. This possibility, 
however, was subsequently withdrawn by the authors due to failings on energetic grounds. In Ref. a high redshift 
{z ^ 10 ) cosmic ray a-burst was postulated and claimed to account for the data but no suggestion about the origins 
for these cosmic rays was made. The sources which may have a cosmic ray fluence sufficient to synthesize such large 
amounts of ®Li were analysed in detail in Ref. Q. It was concluded that typical core collapse supernovae, believed to 
be the source of the standard cosmic rays, fail by a large factor. A similar conclusion was reached for shock-generated 
cosmic rays at the formation of the Milky Way. Only two candidates seemed to fulfil the energetic requirements, with 
both, nevertheless, involving fairly drastic and observationally unconfirmed assumptions. The ®Li may have been due 
to cosmic rays generated by a very early and efficiently forming population of 30-100 Mq stars (involving around 
10% of all baryons in the Milky Way halo), forming black holes at the end of their lifes and ejecting only negligible 
amounts of iron. This scenario implies a large variation of the so far believed universal Salpeter stellar mass function 
and would lead to a large number of supermassive black holes in the galactic halo. Alternatively, the energetic particles 
responsible for the ®Li synthesis may be due to accretion of baryons on the black hole of the galactic center of the 
Milky Way, if the black hole formed before the assembly of the galactic halo and, if accretion on the black hole was 
around a factor of 10^ more efficient in the distant past than observed today [i^. A scenario of this sort, if operative 
in other galaxies as well, may also play a desired role on the entropy in galaxy clusters and contribute significantly to 
the extragalactic 7 background [^. On the other hand, the extragalactic 7 background may have the potential to 
rule out the required ®Li synthesis by cosmic rays altogether, particularly if ®Li has been substantially (factor 10 — 40) 
depleted in PopII stars. Such large ®Li depletion factors are generally predicted when ’^Li is depleted at a factor of 
2.5. Finally, a possible connection between the by WMAP implied fairly early reionisation of the Universe and the 
synthesis of ®Li by cosmic rays has also been considered [4^. 

Well after the submission of the present manuscript to astro-ph, and shortly after publication, it 
-was realized that catalytic reactions invol ving bound states between the electrically charged stau and 
nuclei could considerably change results [^, l50l| . In what follows, to update the paper, all results 
of the original version of the manuscript are shown in the left panels of figures, whereas results fully 
accounting for catalytic effects are shown in the right panels. These are performed by using the recently 
computed reaction rates in Ref. |5ll . l52l| Since the main conclusions of the paper are not changed when 
considering catalytic effects, the text of the remaining manuscript has not been modified from its 
initial version. Prom the figures it is seen that catalytic effects rule out stau decays with life times 
exceeding r « lO^sec due to ®Li overproduction. 
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II. SOLVING LITHIUM PROBLEMS IN THE CMSSM 

It may be more economical to suppose that the ®Li was synthesized in the early Universe, during or right after BBN. 
This may seem particularly attractive if the same process also solves the ^Li discrepancy. In this paper we study 
this possibility in the well-defined Constrained Minimal Supersymmetric Standard Model (CMSSM) coupled 
to gravity, under the assumption that the gravitino is the lightest supersymmetric particle (LSP). In the CMSSM 
supersymmetry is broken in a hidden sector and SUSY breaking is communicated to the visible sector by gravitational 
interactions. The CMSSM is parameterised by five quantities: a unified scalar mass mo, a unified gaugino mass mi/ 2 , 
a unified trilinear coupling Aq (taken to be zero throughout this paper), tan/3 - the ratio between the two Higgs 
vacuum expectation values and sgn(/r) where ji is the supersymmetric mixing parameter of the two Higgs doublets. 
(We assume sgn(/i) = 1.) The first of these quantities are input to the renormalization group equations assuming 
fixed values of tan/3 and are followed from the GUT scale to the electroweak scale in order to determine the 
supersymmetric mass spectrum at the weak scale. Two additional parameters are the gravitino mass mg and the 
reheating temperature Tr, the latter of which is only relevant if TP is efficient. The relic density of the NLSP is 
determined with high accuracy by following the freezeout from chemical equilibrium. We include all annihilation 
and coannihilation channels. For details the reader is referred to Ref. [l^, El. Whereas a relatively fine scan over 
the parameters mo and mi /2 is performed, for the gravitino mass we assume a number of heuristic relations, i.e. 
mg = 0.2mi/2, mg = 0.2mo, mg = mo and mg = 1,10 and 100GeV. 

For each particular point in the parameter space a complete BBN computation is performed using the code as 
introduced in Ref. [l|. This code incorporates all the relevant hadronic and electromagnetic interactions, including 
the most recent data on non-thermal nuclear spallation and fusion reactions, required to make precise abundance 
predictions. A baryonic density of Q.bh? = 0.022 is assumed leading to the following abundances in the absence of 
particle decay: ^Li/Hw 4.31 x 10“^°, D/Hss 2.67 x 10“^, and ^He/Dss 0.39. Apart from the relevant accelerator 
and laboratory limits on particle physics beyond the standard model each point is also subjected to the following 
observationally inferred limits on the light-element abundances: 2.2 x 10“® ^ D/H < 5.3 x 10“® derived from the 
D/H abundance in the local interstellar medium and in high-redshift Lyman-a absorbers, ^He/D ^ 1.72 derived from 
the presolar nebula and, Yp < 0.258, where Yp is the helium mass fraction. 

For the lithium isotopes we examine three conceptually different possibilities: The ®Li is due to relic particle decay 
with the primordial ^Li abundance not much changed. The discrepancy between the standard BBN predicted and 
observationally inferred ^Li/H ratio is solved by stellar depletion. A ^Li destruction factor of 2.5 implies generically [13, 
[Isjl a ®Li destruction of 10 — 40. Therefore we apply: 

(a) ^Li/H > 2.5 X 10-i° and 0.015 < < 3. 

The ^Li/H ratio is considerably reduced by relic particle decay with the ®Li due to other pre-galactic sources. This 
case corresponds approximately to: 

(b) 9 X 10-11 < I'Li/H < 2.5 X 10-1° and °Li/i'Li < 0.015, 

and finally, both lithium problems are solved by relic particle decay corresponding to: 

(c) 9 X 10-11 < '^Li/H < 2.5 X 10-i° and 0.015 < °Li/'^Li < 0.15. 

As the °Li isotope is more fragile than the isotope, it is conceivable that some stellar °Li depletion has occurred 
(e.g. on the pre-main-sequence) even in the absence of '"Li depletion. We have therefore also considered a case 

(d) 9 X 10-11 < I'Li/H < 2.5 x 10-i° and 0.15 < °Li/i'Li < 3. 

where the °Li is in excess of the observations. This case is shown by pink (with grey-shading between the shadings 
of green and red) and is found directly adjacent to the area (c). We note here that we have relaxed the limit applied 
in Ref. [ll| on the acceptable ^Li/i’Li ratio in order to account for significant °Li depletion. 


HI. RESULTS 


Results of our BBN calculations with decaying NLSPs in the CMSSM are shown, for the particular case of tan /3 = 10, 
/i > 0, and Aq = 0, in Figs. 1-4 for different choices of mg. (We assume low enough Tr so that the TP contribution 
is negligible.) Fig. 1 shows the cosmologically interesting parameter space in the mi /2 ~ fno unifying GUT-scale 
mass plane. To help understanding the figures, we remind the reader of some basic mass relations. The mass of the 
gluino is roughly given by rrig ~ 2.7mi/2- The mass of the lightest neutralino, which in the CMSSM is almost a pure 
bino, is ~ 0.4mi/2. The lighter stau fi is dominated by tr and well above mz its mass is (neglecting Yukawa 
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FIG. 1: Parameter space in the GUT-scale unified supersymmetric scalar mass mo - gaugino mass mi /2 plane (all in GeV) 
where NLSP decay into gravitinos may resolve one or both of the lithium problems. The right panel shows results with 
catalytic reactions included, whereas the left panel does neglect such reactions. The parameters of the CMSSM point employed 
are tan /3 = 10, /r > 0, and = 0 and a number of different gravitino mass mg choices as explained in the text. The 
origin of the ®Li in low-metallicity stars may be explained (criterium (a), see text) in the area indicated by green (light grey). 
The discrepancy between observationally inferred- and standard BBN predicted primordial ^Li/H abundance may be resolved 
(criterium (b), see text) in the area shown in red (darker grey). Both lithium problems may be solved at the same time 
(criterium (c), see text) in the area shown by blue (darkest grey). When additional stellar ®Li depletion (see text) occurs, both 
lithium problems may be resolved (criterum (d), see text) in the small area shown by pink (grey shading between the shading 
of green and red). 


contributions at large tan/3) roughly given by cs ttiq + O.lbrriy^- Points which satisfy the lithium abundance 
criteria of case (a), i.e. providing potentially a solution to the origin of the high ®Li abundance in low-metallicity stars, 
are displayed by the green (light grey) points. It is clear that much of the parameter space may solve the ®Li problem 
as this isotope is easily synthesized during a perturbed BBN, either by hadronic decays at earlier times r ^ 10^ sec 
or by hadronic and electromagnetic decays at later times r ^ 3 x 10® sec. If there is any non-thermal and sufficiently 
energetic source during or after BBN, ®Li is normally the first element which is significantly perturbed compared to 
the observations [s^. The parameter space where only the ^Li abundance is significantly reduced, i.e. case (b), is 
shown by red dots (darker grey). This parameter space mostly corresponds to early decay 10^ sec. Finally, both 
lithium problems may be solved at the same time in the area which is shown by blue points (darkest grey). Except for 
a degeneracy in mo this area is well defined, and solutions are found at rni /2 ~ 3TeV. It corresponds exactly to the 
proposed solution to the lithium problems in Ref. [l|, i.e. the decay of a relic of abundance ss 1 — 5 x 10“"^ 

at close to 10® sec after the Big Bang. 

The shape of the area which solves both lithium problems (blue) in Fig. 1 is actually dependent on our discrete 
choices of the gravitino mass parameter. Here the vertical band corresponds to 100 GeV gravitinos, whereas the lower- 
and upper- horizontal bands correspond to the choices mg = mo and mg = 0.2mo, respectively. If we were to vary 
the gravitino mass as a completely free parameter the blue area would be thus significantly enlarged. 

We alert the reader that due to our consideration of multiple possibilities for the gravitino mass, for given mo and 
mi/ 2 , the parameter space in Fig. 1 (as well as in Fig. 2-4, see below), shows simultaneously several choices for the 
gravitino mass. This implies that, for example, points which show that criterium (c) is satisfied may cover up to 
coincide that for the same ijiq - mi /2 but a different mg also a point satisfying, for example, criterium (a) may be 
potentially found. We have plotted, from bottom layer to top layer, first all points satisfying constraints (a), then 
(b), (d), and (c), such that all points satisfying (c) are visible. 

It is encouraging that the CMSSM coupled to gravity provides solutions to both lithium problems. Essentially all 
of these solutions are obtained in the parameter space where the stau is the lightest ordinary supersymmetric partner 
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FIG. 2: Present day gravitino abundance as a function of NLSP decay time in the points shown in Fig. 1. The right 

panel shows results with catalytic reactions included, whereas the left panel does neglect such reactions. The color coding 
is that of Fig. 1. Here only the gravitino abundance generated during NLSP decay (the NTP component) is shown. An 
additional contribution to the gravitino abundance could result for a sufficiently high cosmic reheat temperature T ~ lO^GeV 
after inflation. 



FIG. 3: Gravitino mass mg as a function of NLSP mass mNLSP (all in GeV) for those points shown in Fig. 1 and 2. The right 
panel shows results with catalytic reactions included, whereas the left panel does neglect such reactions. The color coding is 
that of Fig. 1. 
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FIG. 4: Present day free-streaming velocity vo of the gravitino dark matter generated during NLSP decay as a function of the 
fractional contribution of gravitinos to the critical density in those points shown in Fig. 1, 2, and 3. The right panel 

shows results with catalytic reactions included, whereas the left panel does neglect such reactions. Color coding is explained in 
Fig. 1. No reheat-temperature dependent thermal production of gravitinos has been considered. 


particle and the NLSP. In the CMSSM alone, this parameter space is often claimed to be cosmologically disfavored 
due to the electric charge of the stau. This is in stark contrast to our findings that points in the stau NLSP region 
should be regarded as cosmologically favored due to their effect on BBN. Parameter combinations where the bino 
is the NLSP are disfavored, when seen in light of the cosmic lithium problems, as the bino freezeout abundance is 
usually appreciable ~ 1 and it’s hadronic branching ratio is large due to decay into the Z boson, thus providing 

too strong ^ 5 x 10“"^ of a perturbation to BBN. 

Decay times and final present-day gravitino abundance of the points shown in Fig. 1 are shown in Fig. 2. Here only 
the NTP component of gravitinos produced during NLSP decay is shown, adequate for a low reheat temperature Tr 
after inflation. It is seen that the cosmologically most appealing points occur indeed for NLSP decay times around 
lO^sec. Furthermore, it is intriguing to note that those points (blue and pink) at the same time may provide a 
significant fraction, or all, of the by WMAP required dark matter density ~ 0.11 in form of the created 

gravitinos during the decay. This favorable coincidence is related to the fact that the hadronic branching ratio of 
the stau is usually small Bh « 10““^ — 10“^ since a higher-order process than the dominant decay into a tau and 
gravitino. To obtain Q.fh?Bh « 5 x 10“"^ a stau freezeout density of « 0.5 — 5 is required. Given the stau 
decay time Tf w 0.58sec(mg/lGeV)^(TOf/lTeV)“® and a typical stau mass of 1 TeV the desired decay time 10^ 
sec is obtained for a gravitino mass of mg « 50 GeV, i.e. at the electroweak scale. The gravitino relic abundance 
= Dfft,^(mg/mf) 0.025 — 0.25 then comes naturally very close to that required by WMAP. This may also be 
seen in Fig. 3, which shows the stau- and gravitino- masses for the points shown in Fig. 1. The favorite region (blue 
and pink) obtains for stau masses around 1 TeV and gravitino masses between 30 and 200 GeV. Note that in our 
study in Ref. [HI we had not found regions where the dark matter abundance due to only the NTP (decay-produced) 
of gravitinos may account for the totality of the dark matter, since we had constrained our analysis to mi /2 3 TeV. 
Points which only explain the origin of ®Li are obtained for lighter staus (smaller mi/ 2 ), whereas points which only 
solve the ^Li discrepancy are found for heavier staus (larger mi/ 2 ). 

We note here that naively one would have expected to find additional solutions to both lithium problems for decay 
times around r 10® sec and correspondingly lighter staus. These solutions could occur since ^Be, which provides 90% 
of the primordial ^Li since later electron-capturing, could be photodisintegrated in a narrow temperature window 
without the photodisintegration of any other element (in particular D). This would be due to the particular low 
binding energy of ^Be. If this happened, for the accidentally right hadronic branching ratio Bh ^ 5 x 10“® the 
^Be could be reduced, and ®Li could be produced at the observed level by the small fraction of hadronic decays. 
Nevertheless, experimental data on the photodisintegration process ’’Be( 7 , Q!)®He does not exist. If the cross section 








FIG. 5: BBN yields of ’’Li/H and ^H/H for the CMSSM parameter space shown in Fig.1-4. For color coding cf. to Fig. 1. The 
right panel shows results with catalytic reactions included, whereas the left panel does neglect such reactions. 



10' F 


lO-'" SxlO-'o 3x10-'° 4x10-'° 

m/H 


10 ° - 


10 - 


> 10-2 


10-3 ^ 


10 - 


10 - 


1 ' ' ' ' I ' ' ' ' I 
tan/S=10, /U.>0, Ag=0 



_L 


_L 


_I_I_I_ I _1_I_L_ 


_L 


10-'° 2x10-'° 3x10-'° 4x10-'° 


’’U/H 


FIG. 6: BBN yields of ®Li/^Li and ^Li/H for the CMSSM parameter space shown in Fig.1-4. For color coding cf. to Fig. 1. 
The right panel shows results with catalytic reactions included, whereas the left panel does neglect such reactions. 


for the mirror reaction ^Li( 7 , a)^H [HH ] is used solutions may be indeed found However, when reverse reaction rate 
data is used for the well-studied ^He(a, 7 )^Be reaction the ^Be photodisintegration cross section is found abnormally 
low, such that in practice, ^Be photodisintegration is always accompanied by some observationally unacceptable D 
photodisintegration. 

Gravitino dark matter which is generated by the decay of NLSPs is necessarily warmfish), i.e. endowed with free- 
streaming velocities which impact the formation of structure in the Universe [571ISSL l59l IbOL l6l|| . It is known that 
warm- and mixed- dark matter is constrained by a successful small-scale structure formation and a successful early 
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FIG. 7: As Fig. 1, but for tan/3 = 50. The right panel shows results with catalytic reactions included, whereas the left panel 
does neglect such reactions. 


reionisation. Limits from the Lyman-a forest [i^ and the presence of a supermassive black hole at high redshift 
are typically close to Vrms.o ^ 0.10 km/s for the present-day root-mean-square free-streamin g ve locity. Potentially 
even more stringent limits may be derived by the requirement of early cosmic reionisation |64l . l65j | . These may be as 
strong as Vrms,o ^ 0.03 — 0.002 km/s depending on the exact reionisation epoch, i.e. z ~ 17 as indicated by WMAP 
or 2 ss 6 from high-redshift quasar absorption line systems, as well as on the efficiency of star formation. On the 
other hand, warm dark matter has also been claimed to have beneficial effects on structure formation [^ . such as 
the suppression of small-scale structure in Milky-Way type halos and the introduction of constant density cores in 
dwarf spirals. 

In Fig. 4 we show the present-day free-streaming velocities for the gravitino dark matter generated during NLSP 
decay. It is seen that points in the preferred region which come close to the dark matter density inferred from 
cosmological observations have Vrms,o ~ 0.007 km/s. This corresponds to the free-streaming velocity of a thermally 
generated gravitino of mass 3.7 keV, or equivalently, a reduction of the primordial dark matter power spectrum 
compared to that of cold dark matter by a factor two on a scale of 50kpc, for example. At present such velocities 
do not violate any cosmological constraints, nevertheless, are of a magnitude which may be interesting to small- 
scale structure formation. Further information about the cosmic reionisation history could constrain such points. 
Comparatively high free-streaming velocities are reached for points which satisfy criterium (a). Here the dark matter 
abundances generated by NLSP decay are, however, fairly low, such that most of these points are not ruled out. If 
accompanied by a component of cold dark matter (e.g., TP of gravitinos during reheating) such points become mixed 
dark matter. Constraints on mixed dark matter points from reionisation have recently been re-analyzed in Ref. [^ . 
It is interesting that the considered scenarios provide an additional verifiable/refutable prediction on the warmness 
of the dark matter. 

BBN abundance yields of D/H, ^Li/H, and ^Li/®Li for the points shown in Figs.1-4 are shown in Figs. 5 and 6. 
Fig. 5 illustrates that ^Li/H ratios as low as 1.5 x 10“^° may be synthesized (corresponding to a ^Li “depletion” as 
much as 0.46 dex) while producing an observationally satisfying ®Li/^Li ratio as high as 0.04. On the other hand, as 
already noted in Ref. [l| the same scenarios also lead to an increase in the D/H abundance. The predicted higher 
D/H > 3.5 X 10“® fares actually less well with observations [1^ than the prediction of a SBBN scenario at the WMAP 
determined baryon density. However, interpretation of the available data has to be performed with caution, as there is 
actually a dispersion in the inferred D/H ratios in different low-metallicity Lyman-a absorbers which is much larger 
than the inferred errors in individual D/H determinations. This indicates possibly large and unknown systematic 
errors as the naive expectation would be to find D/H constant in different absorbers. Furthermore,D/H is essentially 
always destroyed in stars, leading to the possibility of a D/H underestimate when much gas has been cycled through 
very massive stars. 
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We have also considered CMSSM scenarios at tan/3 = 50, /r > 0, Aq = 0. Results for this choice of tan/3 are shown 
in Fig. 1. It is seen that, though the origin of the ®Li isotope may be explained, there is an absence of points resolving 
both lithium problems simultaneously. That is in contrast to the case tan/3 = 10. 

We have so far concentrated only on gravity-mediated SUSY breaking where the gravitino mass is expected at the 
electroweak scale. When SUSY breaking occurs in a hidden sector which communicates with the visible sector via 
gauge-interactions the gravitino is the LSP and it’s mass may be rather small 1 keV $ mg ^ 1 GeV. To obtain an 
NLSP lifetime of 10^ sec for, e.g., a 100 MeV gravitino, requires a 90 GeV NLSP, close to the experimental bound in 
case of stau NLSPs. Much lighter gravitinos would require NLSPs with mass in conflict with LEP data. This mass 
pattern would lead to much larger free-streaming velocities, i.e. Vrms,o ~ Ikm/s. On the other hand such a case 
would typically only lead to a very small gravitino dark matter density due to the small mg. These points are thus 
not ruled out and may also solve the lithium problems. However they would not provide the bulk of the dark matter. 
A mass scale of m^visp 100 GeV implies squarks and gluinos typically in the several hundreds of GeV up to some 
2-3 TeV range. Supersymmetry at that scale is expected to be discoverable at the LHC. In contrast, a solution of 
the lithium problems in the CMSSM typically implies squarks and gluinos in the several TeV range, unlikely to be 
produced at the LHC. From the point of view of fine tuning, this range is somewhat less attractive and for this reason 
has not been explored in EP. 

Lastly, we mention also that for much of the interesting parameter space in which the stau is the NLSP, the 
scalar supersymmetric potential includes global minima with energy lower than that of the Fermi vacuum [^ . The 
Fermi vacuum would be thus rendered metastable to decay towards the global vacuum. Since it may, however, be 
cosmologically sufficiently long lived, such configurations are not ruled out as long as the Universe settled into the 
proper vacuum after inflation. Alternatively, when Aq ^ 0 the Fermi vacuum is often the true vacuum, and such 
consideration do not play a role. 


IV. SUMMARY 

In conclusion, the observationally inferred primordial ^Li abundance is a factor 2 — 3 lower than that predicted 
by a standard BBN scenario at the baryon density as inferred by WMAP. Though it is conceivable that ^Li, which 
is observed in the atmospheres of low-metallicity stars, has been destroyed in these stars, there exist currently no 
self-consistent and physically motivated scenarios which may explain the observational data. In contrast, the ®Li 
abundance [13 inferred in low-metallicity stars, an isotope which is usually not associated with BBN but rather 
galactic cosmic ray nucleosynthesis, is substantially larger than those predicted by galactic cosmic ray scenarios. It 
has been shown that only under extreme assumptions about putative early cosmic ray populations may the ®Li of such 
a magnitude result [i^. Both, the ^Li and ®Li problems are linked, as significant destruction of ^Li in low-metallicity 
stars (a factor of 2.5) typically implies an even larger destruction of the more fragile ®Li (a factor of 10 — 40), making 
the required synthesis of ®Li by energetic particles even more problematic. We investigate here in the context of the 
GMSSM, and under the assumption of the gravitino being the LSP, whether one or both of the lithium problems may 
be solved by NLSP decays into gravitinos during or after BBN. Here the NLSP density is computed self-consistently 
and with high accuracy. We have found that there exists ample of supersymmetric parameter space where the origin 
of ^Li may be explained by stau decay during or after BBN. This had been already shown earlier in the context 
of hadronic decays [13 and electroma gne tic decays [l^]. Similarly, the ^Li may be effectively reduced when staus 
decay during BBN. It has been shown [l|] that both problems may be solved simultaneously, by the decay of a relic 
particle at about 1000sec after the Big Bang. These exist for tan/3 = 10 with 1 TeV staus decaying into 50 — 200 GeV 
gravitinos. They also lead to an enhancement of the primordial D/H ratio as compared to that predicted in SBBN. 
By chance, these scenarios result in gravitino abundances produced during decay, H^rpp, which contribute a large 
fraction, or all, to the by WMAP inferred dark matter density Hdm- In cases where < Qdm the gravitino 

density may be additionally augmented by production of gravitinos at reheat temperatures of ~ 10®GeV. A prediction 
of these scenarios is the dark matter to be either lukewarm (for H^rpp = VLdm) or mixed (for H^rpp < Qdm) with free- 
streaming lengths of a magnitude interesting to structure formation and relevant to the galactic core- and substructure 
problems. Scenarios of this sort are potentially constrainable/verifiable by the process of reionisation, the Lyman-a 
forest, and/or weak lensing. It would indeed be interesting if the anomalies inferred in the abundances of the lithium 
isotopes would present us with information about the nature of the dark matter and physics beyond the standard 
model. 
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